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Abstract 

This letter is devoted to results on intermediate asymptotics for the heat equation. We study the convergence towards 
a stationary solution in self-similar variables. By assuming the equality of some moments of the initial data and of the 
stationary solution, we get improved convergence rates using entropy / entropy-production methods. We establish the 
equivalence of the exponential decay of the entropies with new, improved functional inequalities in restricted classes 
of functions. This letter is the counterpart in a linear framework of a recent work on fast diffusion equations, see |@]. 
Results extend to the case of a Fokker-Planck equation with a general confining potential. 
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Consider the heat equation in the euclidean space, 

du J 

— = Au t>0, xeW (1) 
dt 

with an initial condition uq e L'(M''). By writing u - u+ - m_ where u+ and «_ are respectively the positive and 
negative parts of u and solving ([T]i with initial data (mo)+ and (mq)-, we may reduce the problem to the case of a 
nonnegative function, corresponding to a nonnegative initial condition mq, without restriction. The heat equation 
being linear, we can assume without loss of generality that mq is a probability measure so that in the sequel of this note 
J^jUodx - 1 - u{t, x) dx for any f > 0. Getting decay rates and even an asymptotic expansion for large values 



of t is completely standard, see for instance u3\\ . However, a few details and some notations will be useful for later 
purpose. 

First of all, as a straightforward consequence of the expression of the Green function, G(f, x,y) :- (4nty'^^ e^^^, 
any solution m of ([TJ can be written as u(t, x) - J^^ uo(y) G(t, x, y) dy and therefore uniformly decays like 0(t^''^^) 
since, as f ^ oo, u(t,x) ~ G{t,x,0). It is also classical to estimate the decay of u{t, ■) - G{t, -,0) in various L''{M.'') 
norms. Such estimates are called intermediate asymptotics estimates. The point is to determine the first term of 
an asymptotic expansion of the solution as f — > oo. For instance, as we shall see below, it can be proved that 
Mt, ■) - G(t, ; 0)\\om = 0(t-"^) as f ^ 

The entropy method can be used among various other approaches to obtain such an estimate. It relies on the 
logarithmic Sobolev inequality and goes as follows. First consider the time-dependent rescaling 



M(f, x) = RT'' V (log R, x/R) with R = Rit) := Vl +2f , f > , x 6 M'' . (2) 
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If M is a solution of ([T]i, then v solves the Fokker-Planck equation 

dv 

-=Av + V-(xy) (3) 
at 

with same initial condition v(f = 0, ■) = "o- Let v^oix) :- (27r)"''^^ e"'^''^^ be the unique stationary solution of Q 
with mass 1, and define dfi :- Vcodx as the Gaussian measure. We denote by L''(M'') and L''(W', djj) the Lebesgue 
spaces corresponding respectively to Lebesgue's measure and to the Gaussian measure. Understanding the interme- 
diate asymptotics for u amounts to study the convergence of v to Voo, as r ^ oo. Define the entropy by £i[w] :- 
^jW\ogwd^. Let V be a solution of Q and define w(f, ■) := v(f, O/^co, wo w{t = 0, •)■ Then ^£i[w(f, ■)] - 
—Ii[w(t,-)] where Ji is the F/i/ier /«/ormaf/on defined by Ji [vv] := w|Vlogw|^c///. Gross' logarithmic Sobolev 
inequality exactly amounts to £i[v/voo] < \ Ji[v/vco] and so, it follows that 
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£i[w(f,-)] <£i[wo]e" Vr>0. 



By the Csiszdr-Kullback inequality, see for instance llTll . we get ||v(f, ■) - Voo||^, ,„rf, ^ i £i [w{t, ■)] and deduce that 



l|v'(f,-)-Voo|lL.(R<') < ^ y/&dwo]e ' Vf>0. 
Undoing the change of variables ^ and observing that u^dt, x) :- Rity' Vk, (x/R{t)) = G(t + 1 /2, ■, 0), we finally get 



iiM(f,-)-«oo(f,-)iiL'(E<')<^ vf>o, 

which establishes the claimed estimate, namely; ||M(f, ■) - G(f, x, 0)||/,i(K<i) < 0(r'^^) as f — » oo. Such an estimate 
is quite classical. The above method is known as the Bakry-Emery method or entropy / entropy-production method 
and also provides a proof of the logarithmic Sobolev inequality. See iflil [st] for some references on this topic, in the 
context of partial differential equations. 

By combining L'(R'') and L°°(M.'') estimates using Holder's inequality, we get that 

||M(r,-)-G(r,-,0)|b(R.)<o(rA<'^(''-'"'') as r ^ cx, . 

In a L^(R^) framework, a much more detailed description can be achieved using a spectral decomposition. If v is a 
solution of Q, then w = v/voo is a solution of the Ornstein-Uhlenbeck equation 

— = Aw - X ■ Vw (4) 
at 

with initial data wq — mq/Vco. Notice that wodjj - 1 and, as a consequence, w(t, ■)djj - 1 for all t > 0. Define 
by (Hk)iceN'i the sequence of Hermite type polynomials (see for instance ifioll ) acting on x = (xi,X2 . . .x^) € W\ 
such that //^.(x) := U^i hkj(xj) where h„(y) := (-1)" (n!)-'/^ er/2^(g-v-/2)^ y 6 R and ;t = {h,...,kd) € N''. 

These functions provide an orthonormal family of eigenfunctions in L^(W' , dfi) which spans the eigenspaces of the 
Ornstein-Uhlenbeck operator, that is - (AHt - x ■ VHk) - \k\ H^, where \k\ :- ^j- Up to a scaling, {hn)n(z^ is the 
usual family of Hermite polynomials on K. 
If Wo satisfies the orthogonality condition 



r 



woHkdfi^O V/teN''suchthatO< <«, (5) 
then an improved rate of convergence follows, in the sense that 

lIlLW.rf/.) ^^""'11^0- l|b(R^,rf;,) Vf >0. 

If ^ initially holds, we indeed have w(t, ■) H^d^ - for any t > and any k e W'such that < < n. Then, 
since j,\\w{t, •) - l||^2(Rd = ~2 j^j |Vw(f, dfi, the conclusion holds using the following result. 



Proposition 1 (Improved Poincare inequality). Assume that w e L^(M.'') is such that J^^wdj^ - 1 and the condition 
jjjj w Hkdju - holds for any k e N'^ such that < < n. Then the following inequality holds, with optimal constant: 

The proof is no more than a straightforward rewriting of the Rayleigh quotient ||Vw||^,^gj d/j)^^^^ ~ ^'^^i?-(W.d ) ^^'^^^ 
the appropriate orthogonality condition. Notice that polynomials Hk are of degree \k\ so that the Condition can be 
rephrased in terms of moment conditions. See Il3lll4ll for further results in this direction. 

It is natural to search for improved estimates of convergence also in U'iW') with p e [1 , 2) by looking for improved 
functional inequalities whenever condition (|5]l is fulfilled. We may for instance quote [|2t] in which improvements on 
the constant, but not on the rates, have been achieved for p - I- 

For any p e (1,2], consider the generalized entropy 



r wP -1 

£p[w] := —dp 

Js.1 p-i 



This definition is consistent with the definition of £i because, under the condition ^^wc/yu = l,£p[w] - j^^-^j^dp 
£i[w] as /? ^ 1. The functional &p controls the convergence in L''{W' , d/j) using a generalized Csiszar-Kullback 
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inequality. In M, it has been proved that ||w - 111^,^^^,^^^ < ^I^'p max{ |H|^/„, l}£p[w], for any € [1,2]. 
Since ||w||z.i(R<f.rf^, = 1, we have 1 < l|w||^ = 1 + (p - l)&p[w], and so 



l|w - l\\D'(R^,d^,) < {&p[w]) with Jlpis) := ^ [l + {p - 1) s]''"'' ^^s . (6) 

Next, assume that wH^dp - for any k e N'' such that < < n and consider the generalized Poincare 
inequalities, with p € [1,2], namely 



Sp[w] < Sn,p \ |Vw''/2 f dp V w e H\W', dp) . 



(7) 



(9) 



Such inequalities have been established for « = 1 by W. Beckner in [51] with optimal constant Si^p - 2/p for the 
Gaussian measure. By the same method, it has been shown in \X] that for a larger class of measures dp, if (|7]i holds 
for p - I and p - 2, for some positive constants S„j and S„,2 respectively, then it also holds for any p e (1,2) with 

S„.p = ^ [l - ((2 - p)/p)^"'/(2s„,) j _ 

By the logarithmic Sobolev inequaUty and the improved Poincare inequaUty, see Proposition[T] we know that < 2 
and S„_2 = l/«. Hence it follows that Sn.p < [1 - ((2 - p)l p)" ] \- On the other hand, as in {stl, if w is a solution 
of (IHi, then 

^£,[w(f, ■)]=-- r \^wPi^fdp. 

at P jRi 

If Q is satisfied, we conclude using (|7]i and ^ that any solution of ^ with initial data wq satisfies 

Sp[w(t,-)]<6p[wo]e-^^''"'P^' and \\w{t, ■) - l\\v'(R^,d,) < ^p{&pM) e-'^"'"'" Vf>0, 
with A(n, p) ;= I n (p - 1) [1 - ((2 - p)/p)"Y^ ■ The last estimate holds because, for any f > 0, 

\\w(t, •) - l|b,(K.,rf^) < Jlp {6p[w(t, ■)]) < Jlp {&p[wo] e-^"-"'P^') < Jlp {6p[wo])e-'^"-P^' 

Notice that A{l,p) = 1 and A{n, 2) - n. Nothing is gained as — > 1, since lim^^i A{n, p) = 1 is independent of n. 

On the other hand, by Holder's inequality, we have for free that ||w - l\\ij>(K',dii) ^ 11*^ ~ M\i?-(«.'',dp)- Hence, if w 
is a solution of with initial data wq, we know that \\w(t, ■) - l||z/(R<',rf;i) < e^"' \\wo - l||L2(R<i_d^) as / ^ oo, for any 
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p e [1, 2], if (|5]l is satisfied. By interpolation, we recover the rates of lfl3[[l4ll . However, this is not satisfactory since 
neither ||wo- l||L/'(K'',rf;j) nor £p[wo] are involved in the right hand side of the above estimate. 

Consider first the case p - I- An alternative approach is suggested by the method of I?! @]j which applies to 
the fast diffusion equation ^ - Am'" for m < 1 . By assuming some uniform bound on the initial data, which is 
preserved along the evolution, it is possible to relate the asymptotic rate for intermediate asymptotics with the spec- 
trum of the linearized operator We can indeed observe that ||wo - l|l2,2(]jrf ^ W^o - mLH»',df2) W^o - l|lL~(K'',i//j) ^ 
I Vfii [wo] ||wo - 1||l~(K'','/a') using Holder's inequality and the Csiszar-Kullback inequality. This proves that 

\\w{t, ■) - lll^nKrf,^^) ^ 5 IK - ^L^^(R'',d^l) ^J&\[m] e'"' as f ^ oo 

if (|5]l is satisfied initially. Still, this provides neither an estimate of J^j w{t, ■) log w(t, ■)dii nor a functional inequality 
which improves upon the logarithmic Sobolev inequality. To prove such an inequality, we keep following the strategy 
of 10]. A simple but key idea is to observe that the functions defined for any p e [1, 2] by hp{Q) - 1, hp{\) - p/2 and, 
for any s e (0, l)U(l,oo)by/!^(s) := [s^ - l-p{s- l)]/[{p-l)\s-l\-] if p > l,hi{s) := [s logs-{s- l)]/\s- Ip, 
are continuous, nonnegative, decreasing on and achieve their maximum at 0. Define on L°°(M'') the functional 

'KpM := llwll^^^'l'jj.j sup/ip(w(x)) = l|w||^;^jj^j /zpQnf . 

Theorem 2 (Improved logarithmic Sobolev inequality). Assume that w e L;j°(R'') is such that J^jwdfx = 1 and 
satisfies the condition w Hi; d^ — Ofor any k e N'^ such that < < n. Then the following inequality holds, with 
optimal constant: 

r , , , 'HiM r ivwp 

( w log wdfj. < I dp . 

Proof. We may indeed observe that by the Poincare inequality and using the definition of "Hi, we get 

riz^flf^>^ — r \vw\^dp>-^ — r \w-i\-dp>-^ r w\ogwdp. 

The optimality of the constant can be checked by a lengthy but elementary computation using the functions :- 
Hk(x)x{xs^^^^"^) + Cg for some smooth truncation function such that Q < x ^ l,Af=lon ^(0' 1) and;^^ s in 
M'' \ B(0, 2). Here for ^ e N'' is such that \k\ - n and the constant is chosen so that j^,, w'^^^dp - I. □ 

As a consequence of the Maximum Principle applied to the heat equation ([T]i and the fact that to uq = Vco corre- 
sponds a self-similar solution of ([TJ, namely u(t, x) = G(t + j,x, 0), we have the estimate 

•KiMf,-)] <"Hi[wo] Vf>0. 

By applying Theorem|2] we obtain a new result of decay for £i [w(t, ■)] with a constant which is exactly S\ [wq], to the 
price of a rate which is less than 2 n. 

Corollary 3 (Improved decay rate of the entropy). Let w be a solution of Q with a nonnegative bounded initial 
data wq e L'(R'', dp) such that w^dp — I and (|5]l is satisfied. Then 

£i[w(f,-)] <Si[wo]e""'^'^''"'°' Vf >0. 

This result is actually equivalent to Theorem[2] as follows by differentiating the above inequality at r = (for which 
equality is trivially satisfied) and using the fact that - jj^,, |VwoP/woii/U = ^fii[w(f, ■)]|f=o < £i[wo] ji e""'^'^'f"'°^|,=o- 
What we have achieved is a global, improved exponential decay of the entropy £i in a restricted class of func- 
tions. To simplify even further, for any e e (0, 1) and n 6 N*, consider the set X'l :- [w e V-{W' , dp) : 
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l-e<vv< 1+s a.e. and the condition J^^wH^di^ - holds for any keN'' such that < \k\ < n }, which is ap- 
propriate to handle the optimality case corresponding to e 0+. The best constant in Theorem |2] is indeed asymp- 
totically equivalent to the sharp rate of convergence in Corollary |3] in the sense that limE^o+ infne,v2 n/'?^i[vv] - 
lim^^o, «/[(! + e) hd - £)] = 2 «. 

For simplicity, we have considered only the case p = I, but the method also applies to any p e (1,2). We 
obtain an improved version of © under the restriction that w € L'(R'', dfi) is bounded nonnegative and the condition 
J^.wHkd/j - holds for any A: € N"' such that < < n. With S„,i - 4'7Yi[w]/« and S„,2 - 1/n, we get 

Sn,p < 'K[n,p, w] :- (n (p - 1))"' [l - ((2 - )^"^''"'^] by (O. Using the entropy / entropy -production identity Q, 
the fact that "7C[n, p, w(f, ■)] < 'K[n, p, wq] and the generalized Csiszar-Kullback inequality we obtain 

£p[w(f,-)] <£p[wo]e"^^^» and ||w - l|b(R</, < (£„[wo]) Vf>0. (10) 

Alternatively, an elementary computation as in the proof of Theorem|2]gives a similar result; 



^ r \vw"'^\-dti^ f w"-^\vw\^d,,> — 1^ r \vw\^dfi> — ^ r \w-ifdfi>-^ 



if w c/ju = 1 and the condition J^^w Hi^diu - holds for any keW' such that < \k\ < n. This proves that 

Using (|9l) and Q, this proves that any solution of Q with initial data in wo € L' n L'"(]R'', c///) satisfies 

£p[w(f, ■)] < £p[wo] e-'"^ and ||w - l||z,,.(R^',<i^) < (£p[wo]) e-«/"/(2«.[-ol) y r > . (11) 

Comparing the rates of ( flOl i and (fTTT i is a natural question. In the limit e — » 0, inf,ve;^j "J^plw] ~ sup„,g y,, 'Hp[w] p/2 
and it follows that limg^o —nrT' — rr = - « (p - 1)/[1 - ((2 - p)l p)''] < 2n - limg^o n^V. ^ - Hence, at least in the 
regime e — > 0, (fTTT i is a better estimate in terms of rates than ( fTOb . Undoing the change of variables (|2]), we have 
achieved a detailed result on improved Mo- 
Corollary 4 (Improved intermediate asymptotics for the heat equation). Let p 6 [1,2] and assume that uq is a 
probability measure such that wq — uq/Voo is bounded and satisfies the condition uq dx - Ofor any 6 N'' such 
that < < n. Ifu is the solution of ([T]l with initial condition uq, then 

\\u{t, ■) - uUt, ■)\\Lr(m < (27r)-^ :Hp {&p[wQ]) (1+2 f)- (1-^) V f > . 

i-i 

The proof reUes on the remark that \\u(t, ■) - u^{t, ■)IIz/(R'') ^ l|Moo(f, OII^^Ir,,) WMt, ■) - IWv'iM.'i.df^) where Uoo{t, •) 
G(t + 1 /2, ■, 0). The conclusion holds using Wu^t, •)IIl»(K'') = {27tR-)-'"~ with R = VT+Tf. 

Up to now, we have considered the simple case of the harmonic potential, V(x) = j |xp. As in HI, the previous 
results can be extended to more general potentials as follows. Consider V e Wj'^'^ Pi W^^iM.'') such that J^^ e^^^'^^dx - 1 , 
and define the probability measure dfi(x) := e^^^ '^^dx in M'^', which generalizes the Gaussian measure. Under the above 
conditions on V, the logarithmic Sobolev inequality holds (resp. ^ for p - I) for some positive constant (resp. for 
Si.i > 0). The Ornstein-Uhlenbeck operator N :- -A + W ■ V is essentially self-adjoint on L^idji), has a non- 



degenerate eigenvalue Aq -0 and a spectral gap Ai > 0. According to lllSL Theorem 2.1], N has a pure point spectrum 



without accumulation points. Since limj^^oo A^ - oo, then by 11151 Theorem XIII. 64], the eigenfunctions of N form a 
complete basis of L^{M.'', dfi). We shall denote the eigenvalues by /l^, k eN, and by the corresponding eigenspaces. 
Theorem|2] adapts without changes. Assume that w e L"(M'') is such that J^^wdp - 1. Then 



I w log wdp < — - — I dp 
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under the orthogonality condition: w e (ULi ^kj , that is wfk dfi - Q for any fi; 6 E^, k - 1, 2,. . .« - 1. Next, 
consider the solution w of the Ornstein-Uhlenbeck equation 

^ = -Nw = Aw- Vy ■ Vw, (12) 
ot 

with initial condition wq e (U^=i ^*:) ^ L°°(IR.'') is such that wq dfi - I. With the same definition as above for &p, 
for any solution of ( fT2l ) with initial data wo, (fTTT i is now replaced by 

Sp[w(t,-)]<&plwo]e-''"P"'"ri"-<'^ and ||w - 1|L,.(k.,^^) < y[4£p[wo]) e-''"'"/P«.[-ol) Vf>0. 

Let us conclude this letter by some comments and open questions. It is standard in entropy / entropy-production 
methods that determining sharp rates of convergence in an evolution equation is equivalent to finding sharp constants in 
functional inequaUties, as we have seen in the case of the heat equation: the rate of convergence in L^(M.'', djj) is given 
by the Poincare inequality, while the rate of convergence in entropy, which controls the L'(]R'', norm, is related 
with the logarithmic Sobolev inequality. This is also true for nonlinear diffusion equations, see for instance 1.121 . In 
this case, a breakthrough came from the observation that uniform norms can also be used, see iflQilTilell. to the price of 
a restricted functional framework. This allows to relate nonlinear quantities of entropy type with spectral properties of 
the linearized problem, in an appropriate functional space and, again, to relate sharp rates with best constants, see |@1. 
As long as nonlinear evolution problems are concerned, only a few invariant quantities are usually available: the mass 
and the position of the center of mass of the solution, for instance. In linear evolution problems, we can impose an 
arbitrary number of orthogonality conditions, which are preserved along the evolution. Improved rates of convergence 
are then expected, even when measured with nonlinear quantities like the entropy. Various attempts have been done, 
see for instance 12|], but the question has been left open for many years. Such ideas have been partially explored by 



R.J. McCann, including in the linear case (see |1 based on considerations on an appropriate Hessian matrix. Our 
approach provides a simpler and elementary answer under restrictions which are natural in view of |@] . It also raises a 
number of questions concerning the optimality of the new functional inequalities from a variational point of view, the 
convergence of minimizing sequences and the symmetry of the eventual minimizers. 
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